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j  Abstract 

\  This  paper  presents  approximations  for  the  rms  error  of  the  maximum 
likelihood  estimator  of  the  direction  of  a  plane  wave  Incident  on  a 
random  array.  The  sensor  locations  are  assumed  to  be  realizations  of 
independent.  Identically  distributed  random  vectors.  The  second  part  of 
the  paper  presents  an  asymptotically  unbiased  estimator  of  the  noise 


wavenumber  spectrum  from  random  array  data.  J, 
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Introduction 


Sonobuoy  fields  are  used  to  detect  submarines.  Thorn,  Booth,  and 
Lockwood^  have  proposed  that  the  signals  from  randomly  deployed  sonobuoys 
be  coherently  combined  to  make  acoustic  measurements.  They  present  the 
expected  value  and  variance  of  the  pattern  function,  and  the  distribution 
of  the  directivity  Index  of  a  three-dimensional  random  array.  In  their 
model,  the  sensor  locations  are  observed  realizations  of  random  variables 
that  may  be  correlated  and  have  different  distributions.  They  define  an 
array  to  be  totally  random  If  the  sensor  locations  are  realizations  of 
Independent,  Identically  distributed  random  variables.  Several 
stochastic  properties  of  the  sldelobe  pattern  of  a  totally  random  array 
are  given  by  Steinberg.^ 

The  ratio  of  the  peak  sldelobe  to  the  main  lobe  and  the  directivity 
Index  of  an  array  system  are  measures  of  its  ability  to  perform  Its 
tasks.  The  generic  signal  processing  Casks  of  an  array  system  are: 

1)  detecting  and  estimating  parameters  of  coherent  wave  signals  that 
impinge  on  the  array;  2)  resolving  multiple  wave  signals;  3)  estimating 
range,  bearing,  or  velocity  of  a  source  chat  generates  the  detected 
signal;  and  4)  estimating  Che  frequency-wave number  spectrum  of  Che 
ambient  noise  field.  This  description  of  system  tasks  emphasizes  the 
statistical  nature  of  the  problem  of  measuring  performance,  especially 
for  random  arrays. 
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This  paper  presents  approximations  for  the  mean-square  error  of  the 
maximum  likelihood  estimator  of  the  bearing  of  a  plane  wave  Impinging  on 
a  random  array  from  a  distant  source.  The  second  part  deals  with 
estimating  the  ambient  noise's  wavenumber  spectrum. 

1.  Random  Planar  Arrays 

Consider  a  planar  array  of  M  sensors  where  the  sensor  locations 
{(x^,y)^)}  are  realizations  of  Independent,  Identically  distributed  random 
variables  {(Xj^tY^)}.  Assume  for  simplicity  that  the  signal  Is  a  single 
frequency  plane  wave  plus  stationary,  zero  mean,  Gaussian  noise.  Let  6q 
denote  the  wave's  direction  of  arrival  with  respect  to  the  x  axis.  This 
angle  Is  the  source  bearing  If  the  medium  Is  horizontally  homogeneous. 

Let  Ug,  Xgi  and  A  denote  the  wave's  frequency,  wavelength,  and  complex 
aiqplltude,  respectively.  The  signal  at  the  kth  sensor  is 

s(t,xk,yk)  -  Aexp[i(wQt-K^j^-iCyyj^)]  +  e(t ,xi^,y^) ,  (1) 

where  Xjj  *■  (2ir/AQ)cos0Q  and  Ky  «  (2x/XQ)8in0Q  are  the  x  and  y  components 
of  Is  the  wavenumber,  and  e(t,X}^,y^}  is  a  realization  of  the  noise 
field. 

The  correspondence  between  beamforming  and  frequency-wavenumber 

processing,  and  an  approximation  to  the  maximum  likelihood  (ML)  estimator 

of  0o  have  been  presented  in  a  previous  paper.-'  if  and 

M  -2 

pZk-l'yk"y^^  are  large,  where  p  is  the  power  slgnal-to-nolse  ratio  In  a 

—  1  ^  A 

narrow  band  about  Uq  and  x  ■  ”  Levln^  shows  that  the  root 


mean-*8quare  errors  of  the  ML  estimators  of  and  ic^  are  approximately 


.  M  _  -1/2 

rmse  Xx  “  l^p  J  (xjj-x)^] 
k-1 


M  _  -1/2 

rmse  »  [2p  (yfe-y)^] 

k-1 


(2) 


“  M  _  _ 

Moreover,  the  covariance  la  E(kx  -  Xj^Xxy  -  <y)  *  (2pSyj^*^(xj^-x) (y^^^-y) )” 

These  expected  values  are  conditional  on  a  realized  array  geometry,  l.e 

they  are  ex-post  the  deployment  of  the  array. 

To  approximate  these  errors,  assume  that  M  Is  large.  Since  the 

sensors  must  lie  In  some  closed  and  bounded  set,  the  random  variables 

(X^,Yj^)  are  bounded.  Thus  the  central  limit  theorem  Implies  that 

M  M 

M“lE^.l(Xk-x)2  -  0x2  +  Op(M-^'^)  and  M“lEk.i(yk-y)2  -  ay2  +  Op(M-l/2), 

where  0,^2  and  Oy2  are  the  variances  of  and  Yj^  respectively,  and 

Op(M"1^2)  means  that  for  any  e>0,  there  Is  a  Be>0  such  that  the  error  1 

bounded  by  B£M“l/2  with  probability  1-e.  Thus  the  rms  errors 

m  A 

of  Kx  approximately 


rmse  Xx  *  (2pM)”l/2o“l 

X 

(3) 

A 

rmse  Xy  <*  (2pM)“l /2o"’l 

'  y 


for  large  M.  The  estimators  are  approximately  uncorrelated  If  the 
coordinate  system  Is  rotated  to  make  the  covariance  -  0  after 

xy 

rotation. 


t* 


tt. 
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The  maximum  likelihood  estimator  of  the  bearing  Is 

AAA  A  A 

Bq  ■  tan”^Ky/Kjj)  radians.  The  linear  approximation  of  tan~^(Ky/K^)  - 
tan”^iCy/iCx) 


-2  -1 
(l+<y^Kv  )~^[Kv 


“2 


y  y' 


y'^x 


(Kx‘*''x^  ]  • 


(4) 


Since  and  x  are  approximately  uncorrelated  If  o^  =  0,  It  follows 

A  y  xy 

from  (3)  and  (4)  that  when  pMo^^  and  pMOy^  are  large, 

A 

E(0o“flo)^  a  (Xq/Zti )2(2pM)“’^(o~^sln^6o  o^^cos^Bq)*  (5) 

X  y 

Thus  If  Ojj  *■  Oy  =  a,  then  from  (5) 

^  1/2 

rmaeBo  “  Xq(2pM)“  (2Tro)’'^rads,  (6) 

For  example,  let  o/Xq  ■  12,  M  -  90,  and  p  =  l/4(-6  dB).  Then  from  (6), 
rmseBQ  ■  0.11*  (1.98  x  10“^  rads).  If  o/Xq  “  100,  M  ••  40,  and 

A 

p  ■  -10  dB,  then  rmseBo  ■  0.03*. 

Now  suppose  that  and  are  Independent  uniform  variates  whose 
range  Is  (0,L),  l.e.  the  sensors  are  uniformly  distributed  on  the  square 

A 

{0  <  X  <  L,  0  <  y  <  L} .  Then  o^  ■  L^/12.  Let  us  compare  the  rmseOo  of 

O 

this  random  array  with  that  of  the  square  lattice  array  whose  M-N^ 
sensors  are  at  the  points  {  (Jd,£d)  :  j  ,il-l , . . .  ,N}  .  If  the  length  of  the 
square's  sides  Is  L,  then  the  sensor  spacing  Is  d  •>  L/(N-1). 

From  (2),  (4),  and  (5),  we  only  have  to  compare  M”^Z(xi^-x)^  * 

M"^i:(yk"y)^  with  0^. 


Since 
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M-1  I  (xk-^2 

k-1 


M-^d^N  I  (J-J)2 

J-1 


(7) 


,  l2  i«-1 

-  d'‘(N-l)(I«-l)/12  - - 

12  N-1 

-  l2/12  -  o2, 

expression  (6)  holds  for  the  square  lattice  array.  The  approximate 
rmse  of  the  maximum  likelihood  bearing  estimator  for  a  uniform 

A 

random  array  on  a  square  Is  equal  to  the  approximate  rmseSQ  for  a 
uniformly  spaced  lattice  array  on  the  same  square. 

2.  Three-Dimensional  Random  Arrays 

For  a  given  coordinate  system,  let  denote  the 

vector  location  of  the  kth  sensor  In  a  three-dimensional  array.  Let  6^ 
denote  the  azimuth  angle  of  propagation  with  respect  to  the  x  axis,  and 
let  Oq  denote  the  elevation  angle  with  respect  to  the  z  axis.  Thus  the 
signal  at  the  kth  sensor  Is 

s(t,xjj)  -  AexplKUgt  -£'xn)l  +  e(t,2^), 

where  jc'  ■  (iCx*'^y*'^z^  vector  of  wavenumber  components 

Kj^  ■  (2x/XQ)cos0g,  Ky  -  (2iT/XQ)sln0Q,  and  ■  (2x7x^)0000^. 

The  correspondence  between  beamforming  and  frequency-wavenumber 
processing  holds  In  three  dimensions.  The  ML  estimators  of  the 
wavenumber  components  are  the  <y,  and  tt^at  maximize 


T 
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N  M 

1  I  I 


J-1  k-1 


-  “o’^J  ^  ^ 


2 


(8) 


where  N  is  Che  number  of  simultaneous  discrete-time  observations  of  the  M 

A  A 

channels.^  The  rms  errors  of  ic^  are  approximated  by  (2),  and 

rm8eK2  “  ,  Once  again,  the  ML  estimator  of  the  source 

A  A  A 

bearing  Is  0q  •  tan~^(tCy/iCjj) ,  and  thus  (5)  holds  for  a  totally  random 
three-dimensional  array  of  M  sensors. 

3.  Estimating  the  Wavenumber  Spectrum 

Consider  Che  problem  of  estimating  the  frequency-wavenumber  spectrum 
of  the  ambient,  zero  mean,  Gaussian  noise  field  around  a  random  array. 
Since  an  n-dlmenslonal  array  is  not  much  harder  to  analyze  than  a  linear 
array,  let  ^  , . . .  ’  denote  the  vector  position  of  the  kth 

sensor  with  respect  to  a  fixed  coordinate  system.  Assume  that  the  are 

realizations  of  Independent  random  vectors  {5^  =  > • • • »^kn^ that 

have  a  common  continuous  multivariate  density  f (x) .  Rotate  the 
coordinate  system  so  that  the  covariance  matrix  of  is  diagonal,  and 
for  simplicity  let  ■  (,^2  ^,2^  i.e.  is  the  variance  of  each 

X^jl^  after  rotation. 

Let  e(c,]c)  be  the  noise  at  point  x  at  time  t.  If  the  noise  field  Is 
stationary  in  t  and  jc,  the  covariance  function  Cj.(t,^)  =•  Ee(t+T,x42^)€(t,3^) 
Is  Independent  of  t  and  x.  The  frequency-wavenumber  spectrum  is  defined 
as 


SgCu.O  -  /cg(T,2.)exp[lU’2.  “  “T)ld^, 


(9) 


i 
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assuming  that  Is  absolutely  Integrable.  The  power  spectrum  of  the 
noise  is  S^CUfO). 

Assuming  that  the  channels  are  sampled  at  times  tj^jA  for 

j~0 . N**l,  define  the  discrete  Fourier  transform 

N-1 

{e(^)  ■  5  h  »  If  Sg((»),0)  Is  bandllmlted 

at  n/A,  then  N“^|£(xjt)p  »  A“^Se('*>»0)  for  large  N.^  Let  us  work  with  the 
e(xp)  to  obtain  an  estimator  of  S^(i^,o>)  for  a  given  u,  which  will  be 
denoted  S^C^)  to  simplify  notation.  The  properties  of  the  estimator 
depend  on  the  following  theorem. 


Theorem.  Define  the  n-^llmenslonal  Fourier  transform,^ 

U(£)  ■  lit*!  e(xit)exp(l£'3^) .  Assume  that  D(o)  -  /f^(x)dx  *■  0(o“")  and 
when  je  *  0,  l^(je)|  <  co“”  for  some  constant  c,  where  ■  Eexp(ijc'J^)  is 

the  characteristic  function  of  These  assumptions  hold  for  the 
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11m  (DM^)”^(2it  )”"/ |R(ic)  1^  die  -  D“^(2Tt)“*'/ |())(ic)  |^d£ 

M+o*  “  ~ 

=■  D“l/f2(x)dx  -  1.  (11) 

From  (10),  limM^.»  (DM^)"!  |r(0)  (2  =  D“ll(j)(0)|2  -  D“1  =  0(o").  Thus  (11) 
implies  that  as  M  and  o+<»,  (Dm2)“^  |  R(j^)  1 2+5(j^)  ^  a  Dirac  delta  function.  If 

iil*l2* 

(Dm2)-1r(£i)R*(£2)  -  0(o““)  +0p(M-l/2).  (12) 

These  limit  results  are  used  as  follows: 

M  M 

E[U(£i)U*(£2)]  “I  I  "Jl2'2Sk^^ 

j=l  k=l 

M  M 

»  (2ir)""  I  I  /Se(v)exp[-iv(xj  -  jc^)  ]exp[i(£i'}^j  - 
j-1  k=«l 

from  the  Inverse  of  (9).  Gathering  terms, 

E[U(£i)U*(£2)]  *  (2n)~^/R(Kj  “  v)R*(ic2  ”  v^)Sj.(v^)dv.  (14) 

Thus  from  the  above  limits  and  (14),  limM^(j+oo(DM^)~^E|U(jc)  = 

(  2x )  6  (£-v  )  Sg  ( V )  d V  -  Sg  ( i< ) . 

If  £1  *  £2,  then  llmn^g^„(DM^)“^E[U(£i)U*(K2) 1  =  0  from  (12).  Thus 
U(£i)  and  U(jC2)  are  asymptotically  uncorrelated.  For  finite  M<<a2”,  the 
correlation  is  0p(M”1^2)^ 

This  theorem  provides  a  basis  for  estimating  Sj.(_k).  One  method  is  to 
divide  the  (time)  sample  into  J  segments  of  successive  observations, 

Nj  ■  N/J,  and  compute  U(£)  for  each  segment.  These  Uj(jc)'s  will  be 
approximately  uncorrelated  if  Nj  is  large.  Thus  from  the  theorem. 
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Uj(0  have  a  complex  Gaussian  distribution  for  each  j  (the  noise  Is 
Gaussian),  2(DM^)“^  |  Uj  (<_)  P/Se(£.)  Is  approximately  chl-squared  with  two 

A 

degrees  of  freedom  and  thus  the  variance  of  S^Cjc)  Is  approximately 

J"^s2(k). 
e  “ 

4.  A  Planar  Array  Example 

Continuing  with  the  vector  notation,  suppose  that  the  sensors  are 
uniformly  distributed  on  the  square  {-L/2  <  xj  <  1./2,  -  L/2  <  X2  <  L/2}  . 
Thus  f(^)  ■  1/L^  for  jt  In  the  square,  a\^  *  02^  =•  *  L^/12,  and 

D  *  /f^(x^)dx  “  L“^.  The  assumptions  for  the  theorem  hold  since  D  =  0(a”^) 
and  ••  4(xjX2L^)”^sln(KiL/2)aln(K2L/2)  =  0(a“^) .  Thus 

(L/M)^E|u(tc)  1^  a  Se(i<)  for  large  M  and  L  In  this  example.  The  estimator 
of  Se(£.)  is  then  (L/M)2j-1j:j^  J  Uj(£)p  using  the  time  segmentation  method. 


*Thls  work  was  supported  by  the  Office  of  Naval  Research  (Statistics  and 
Probability  Program)  under  contract. 
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